Timescales of tunneling decay of a localized state 
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Motivated by recent time domain experiments on ultrafast atom ionization, we analyze the tran- 
sients and timescales that characterize, besides the relatively long lifetime, the decay by tunneling 
of a localized state. While the tunneling starts immediately, some time is required for the outgoing 
flux to develop. This short-term behavior depends strongly on the initial state. For the initial 
state tightly localized so that the initial transients are dominated by over-the-barrier motion, the 
timescale for the flux propagation through the barrier is close to the Biittiker-Landauer traversal 
time. Then a quasistationary, slowly decay process follows, which sets ideal conditions for observ- 
ing diffraction in time at longer times and distances. To define operationally a tunnelling time at 
the barrier edge, we extrapolate backwards the propagation of the wave packet escaped from the 
potential. This extrapolated time is considerably longer than the timescale of the flux and density 
buildup at the barrier edge. 

PACS numbers: 03.65. Nk, 03.65.Xp, 03.75.-b 



I. INTRODUCTION 

Tunneling, being one of the most fundamental con- 
cepts in quantum mechanics, remains a source of strong 
theoretical and experimental controversies on the rele- 
vant timescales of the process [H, Q • A typical tunneling 
setting is the scattering process, where a wavepacket is 
reflected from and transmitted through a barrier higher 
than the incident particle energy. In this case, the group 
delay @], defined in terms of the energy derivatives of the 
reflection or transmission phase shift, describes the mo- 
tion of the wave packet peak. It was also found that the 
group delay for particles tunneling through an opaque 
barrier is independent of the barrier width (the "Hart- 
man effect" [|). 

The measurement of the time spent by a tunneling 
particle in the classically forbidden region can be based 
on the approach of Baz and Rybachenko 0, which 
uses the Larmor precession of a particle with a magnetic 
moment in a weak magnetic field or in effective fields in 
solids due to the spin-orbit coupling 0, Q . In classically 
forbidden regions there is not only precession but also a 
rotation Q of the moment into the direction of the field. 
Biittiker and Landauer [To| analyzed tunneling through 
an oscillating rectangular barrier and found an interac- 
tion time which is closely related to the rotation of the 
magnetic moment in a magnetic field 0, [T^] ■ A tunneling 
time has been measured for electromagnetic waves pass- 
ing through inhomogeneous optical structures [llf, [l2[ 
and waveguides [l3ul4j . Interestingly, in graphene, a sin- 
gle layer of carbon atoms, which provides another exam- 
ple of massless particles, the transport occurs via evanes- 
cent waves, and the Wigner-Smith delay is linear in the 
tunneling distance [l5j . 

Brouard, Sala, and Muga, using scattering theory pro- 
jectors for to-be-transmitted/reflected wavepacket parts 
and for localizing the particle at the barrier, set a formal 
framework that unified many of the existing proposals 



pointing out that the multiple time scales correspond to 
different quantizations of the classical transmission time, 
due to the noncommuting observables and possible order- 
ings involved fliij . This clarified the meaning of different 
partitions of the dwell time into transmitted, reflected 
and interference components. Another research line has 
been the investigation of characteristic times for the tran- 
sient dynamics of the wave function (e.g. the forerunners) 
in a plethora of potential configurations and initial con- 
ditions using asymptotic methods (TTj-Qij- For specially 
prepared states, in particular for confined or semicon- 
fined initial waves with a flat density, these transients 
show diffraction in time [H|,[2(|, i.e. temporal oscillations 
reminiscent of spatial Frcsncl-diffraction by a sharp edge 
2C|. Also, the analysis of the partial density of states 
2l| provides another approach to the understanding of 
the tunneling process. 

An important tunneling-dependent phenomenon is the 
decay of a metastable system [22l-|25|. related, e.g., to 
the state ionization in optics and to the discharging of 
a capacitor in mesoscopic physics (26l - [3^ |. Compared 
to the full scattering problem, the decay configuration, 
or "half-scattering" problem, has been frequently consid- 
ered unproblcmatic because of the absence of the trans- 
mitted/reflected wavepacket splitting and the existence 
of a well known characterization in terms of resonance 
lifetimes. In fact, one may still pose classically minded 
questions on the tunnelling time similar to the ones in 
the scattering configuration, however, without obvious 
answers. A particle may wander in the trapping well for 
a while and then escape through the barrier. For such a 
history the lifetime would be a waiting time in the well 
plus a tunneling time. Can this quantity be defined and 
measured in a sensible way? The understanding can be 
based on the analysis of quantum interference among the 
Feynman paths [33[ or on the consideration of the quan- 
tities defined by operational procedures, as presented in 
this paper. 
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In recent years, the techniques of atom ionization by a 
strong laser field and attosecond probing of electron dy- 
namics opened a new venue for experimental studies of 
the tunneling times. The measurement of ionization of 
He atoms [34| holds the promise to observe the tunnel- 
ing delay of electrons in real time. In the experiment 
access to the dynamics at the tunneling time scale is 
gained through extrapolation of long time measurement 
to shorter times by assuming that a particle that has 
escaped through an energetically forbidden into a classi- 
cal allowed region, follows classical dynamical laws. Us- 
ing these laws and scattering data the moment of es- 
cape from the tunneling region into the classical region 
is determined. A related experiment [35[ measures the 
perpendicular distribution of the nascent quantum me- 
chanical momentum uncertainty of the initial state as it 
is revealed by tunneling. 

The aim of this paper is to study the dynamics of the 
tunneling-induced decay of a localized state and to inves- 
tigate to what extent we can extract information on the 
short time behavior from the long time dynamics. We in- 
troduce a simple quantum mechanical model, which sim- 
ulates the experimental measurement of tunneling time 
and allows for opening and closing the barrier |34| . We 
calculate the probability density and flux to study the 
short-term dynamics near the barrier and the long-term 
dynamics far away from it. For an opaque barrier we 
observe a relatively short transient time scale, where the 
tunneling is developed into a quasi-steady decay process. 
The details of the initial state have an important im- 
pact on the transients. At large distances, the flux and 
the density start to grow at long times, reach a peak 
and then decrease by diffraction-in-time oscillations. We 
found that the tunneling time obtained by extrapolat- 
ing the particle motion from the position of the remote 
detector to the right edge of the barrier is not directly 
related to the timcscalc of the outgoing flux buildup. 

II. GENERAL DESCRIPTION OF TUNNELING 
AND MODEL POTENTIAL 



wave function can be expressed as 

f°° ( ik 2 t\ 
< t < t ) = J G{k)cp k {x) cxp f dk(2) 

with 

G(k)= / ^{x,0)c/> k {x)dx, (3) 
Jo 

where 4>k(x) are delocalized real states in the potential 
U%{x) with the energy E = fc 2 /2. We use here fi = in = 1 
units, where m is the electron mass. 
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To study the tunneling dynamics, we consider the po- 
tential U(x,t), infinite at x < 0. At t < the poten- 
tial holds bound states with wavefunction tpj (x) and en- 
ergy Ej] it changes at t = to allow the tunneling, and 
changes back at t = to to its initial form. The entire time 
dependence is 



U(x,t) = 



Ui(x) (t < and t > to), 
U 2 {x) (0<t< t ). 



(1) 



The initial state ^(x,t = 0) prepared at t = begins 
to evolve at t > and the probability to find the elec- 
tron inside the potential decreases. This is similar to the 
ionization of an atom by a strong laser field, which lets 
a valence electron to tunnel through the barrier. At the 
closing time to, the potential becomes Ui(x) again and 
the decay terminates. In the time interval < t < to the 



FIG. 1. The time-dependent potential Ui(x) at t = and 
t > to and U2(x) within the interval < t < to. (a) shows 
that Ui (x) is a step in the positive-half plane, which is from 
x — to x = a± and Uo from x = Oi to infinity, while (b) 
presents a barrier which extends from x — ai to x — a,2 and is 
of height of Uo. The potential is always infinite in the negative 
half plane. 



After closing the barrier at t — to, the potential gets 
its original form, and the wave function can be expressed 
as 

9(x, t>t )=Y / Bjipj(x) cxp [-iEj{t - to)] (4) 

3 

+ / S(fc)^ fc (x)exp - ?fc2(t T to) , 
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where summation is extended over the bound states of 
the initial potential, and 



U 2 (x) 



B, 



^(x, to)<fij(x)dx, 



B(k) = / y{x,t Q )<p k (x)dx, 



(5) 
(6) 



with ifk{x) being the continuum states in the potential 
Ui(x). With Eqs. © and (g]) wc have fully specified the 
time evolution of the escape problem. 
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FIG. 2. (Color online) The matrix element G(k) of the 
initial state with the delocalized states of the potential U2- 
k is the wave vector of the delocalized states. The solid line 
corresponds to the initial state ^(x,Q) — ipo(x), which is a 
bound state of Ux- For this case the barrier in U2 extends 
from 01 = 1 to 02 = 1.4. The dashed line corresponds to a 
sinusoidal initial state 9(x, 0) = V2sin(7ra;) and to a barrier 
in U2 extending from 01 = 1 to 02 = 1.42. In both cases, 
nd = 2 and Uo = 16. 



To specify the model, we assume that at t < the 
potential Ux(x) is a step in the positive half-axis, being 
zero at < x < a\ and Uq at x > 01, as shown in 
FigUJa). At t = the potential changes to U 2 (x), which 
is a rectangular barrier of the height Uq extended from 
ai to a 2 , as shown in Fig. [Ub). The rectangular bar- 
rier being the simplest example of an everywhere finite 
potential, where the wavefunction and its derivative are 
continuous, allows for the exact analysis of the dynamics. 
We put here ax = 1 and measure the time, energy, and 
momentum in the corresponding units. 

The delocalized eigcnstates of the Hamiltonian corre- 
sponding to Ui(x) have the form: 



fk{x) = 



Ci(fc) sin(fcc), 
^2/nsm[qx + 



(0 < x < ax) 

(x > dl) 



(7) 



where q = sjk 2 - 2U , k > y/2Uo~ with Cx(k) and 9x(k) 
determined by the boundary conditions at x = ax, and 
the norm is determined by (ipk^fk) = S(q — q'). 

The eigenstate of the Hamiltonian corresponding to 



C(k) sm(kx), 
fc (x) = { D{k)e-^ x + F{k)e KkX , 
/2j^sm[kx + 9(h)], 



(0 < x < a{) 
(ax < x < o 2 ) (8) 
(x > a 2 ) 



with the normalization condition (<pk'\<Pk) = S(k — k'). 
The coefficients C(k), D(k), F(k) and the phase 9(k) 
satisfy the boundary conditions of the potential U2(x), 
and Kk = \/2Uo — k 2 . In the tunneling regime k < \/2Uo, 
while in the propagating regime k > \/2Uq, and ink is 
substituted by q defined below Eq. (J7J. 



III. TUNNELING DYNAMICS AT SHORT AND 
LONG TIMES. 

A. Decay of the state 

The decay rate is described by the probability to find 
the particle in the well from x = to x = ax defined as 



wx(t) = / V*(x,t)y(x,t)dx. 



(9) 



The lifetime ti, which also can be seen as the dwell time, 
an important time scale introduced to characterize the 
decay rate, is the time that it takes for the relative prob- 
ability in the well Wx (t)/wx (0) to decay to 1/e. For an 
opaque barrier the decay time can be written approxi- 
mately as ti « ATk cxp(2nd), with n = ^2(Uv — Eq) 
and d = a2 — ax, where fco = V2Eq, Tk = 2/ko is the pe- 
riod of motion for a particle in the well, exp(— 2nd) <C 1 
is the tunneling rate, and A is a barrier-dependent coef- 
ficient of order 1. Therefore, nd is a significant physical 
parameter with regard to the decay process and the tun- 
neling time. Analysis of different examples where this 
simple expression for the lifetime does hold can be found 
in Ref. |25j|. 

To study the dynamics, we first set the initial state 
as the ground state of the Hamiltonian corresponding to 
V~x(x) for a typical potential. For example, for Uq = 16 
there are two bound eigenstates: the ground (fo(x) (en- 
ergy Eq) and the excited state <px(x) (energy Ex). The 
energy Eq for &(x, 0) = <po(x) is 3.52. By choosing dif- 
ferent 02, one can modify the transparency of the barrier. 
For example, if nd = 2, corresponding to a 2 = 1.4, this 
barrier is a moderately opaque one. The correspond- 
ing G(k) is shown in Fig. [5J Under the conditions of 
^(x, 0) = <po(x) the lifetime in this potential is t; = 17.5, 
as illustrated in Fig. G3a)- Calculated real and imag- 
inary parts of ty(x,t) after some short transients show 
fast oscillations with the envelope function exp (—t/2ti). 
This behavior implies that in terms of the poles in the 
complex energy plane, the states we consider correspond 
to the simple Breit-Wigner resonances. Detailed analy- 
sis of various types of resonances and their relation to 
different timescales can be found in Refs. (3(| H3|- To 
clarify the influence of the initial state on tunneling, we 
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FIG. 3. (Color online) The probability to find the particle 
in the well from to oi. (a) The solid and the dashed lines 
are for two initial states ^(a^O) = (po(x), one for Co = 16 
and the other for Uo = 10. (b) The solid and the dashed lines 
represent the barriers Uo = 16, and Uo = 10 with the same 
initial state 9(x, 0) = y/2s'm(-Kx). The insets exhibit ln[p(i)] 
at small time scale, where g(t) = wi(t)/wi(0). For all curves 
Kd = 2. 
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FIG. 4. (Color online) (a) Flux at the left edge ai. The 
solid, dashed and dot-dashed lines correspond to the barriers 
Uo = 24, Uo = 16, and Uo = 10, respectively, (b) Flux at the 
right edge a.2- The inset presents the density at ai (solid) and 
multiplied by fifteen at a-z (dashed), with the barrier Uo = 16. 
All lines correspond to the same parameters as those in (a). 
For both plots, nd = 2 and ty(x,0) = ipo(x). 



alter it into the ground state of an infinite- wall potential, 
^(x, 0) = y/2sm(nx). Therefore, E is tt 2 /2 and a 2 be- 
comes 1.42 to keep nd unchanged. As shown in Fig. [2] 
the coefficient G(k) in Eq. §3§ has two contributions, re- 
lated to the resonances corresponding to the bound states 
of the initial potential. For $(x, 0) = <fo(x), the second 
one, corresponding to the first excited state with a fast 
decay, is extremely weak. The presence of more than one 
bound state in the initial potential combined with the 
condition of low transparency of the barrier leads to im- 
portant consequences for the short-time scale tunneling 
dynamics. 

Although all the probabilities decrease exponentially 
at long times, showing the general feature of the pro- 
cess, the ones with &(x, 0) = V2sin(7ra;) in Fig. [3jb) 
oscillate fast, while those with ^(x, 0) = <po(x) in (a) de- 
cay smoothly. This is because $(x, 0) = \/2sin(7rx) has 
a significant contribution of the second anomaly shown 
in Fig. pleading to interference with the "ground state 



B. Short-term dynamics 

In this subsection we address the beginning of the tun- 
neling by studying the short-time dynamics of the flux 



2,1 



4/ 'd^(x,t) dm*(x,t) , ' 



dx dx 
and the density, 

p(x,t) = ^*(x,t)^f(x,t). 
These two quantities satisfy the continuity equation 
dp(x,t) t dJ(x,t) ^ 



(10) 
(11) 



dt dx 
Based on Eq. (|T2"j) . we obtain the flux at the edges, 



(12) 



J(a lt t) 



dwi(t) dw 2 (t) 
J{a 2 ,t) - 



dt 



where 



resonance 



w 2 {t) 



dt 



^*{x,t)^{x,t)dx 



(13) 
(14) 
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FIG. 5. (Color online) (a) Flux at the right edge on time with 
different width: solid line, nd = 1.5, dashed line, nd = 2, and 
dot-dashed line, nd — 2.5, Uo = 16. (b) Flux at the right 
edge on time with different height, by remaining ai = 1.4: 
solid line, nd = 1.5, dashed line nd — 2, and dot-dashed 
line, hid — 2.5, respectively. The ground state energy Eo 
depends only weakly on the potential height Uo. For all plots 
*(x,0) = <p (x). 



is the probability to find the particle outside the poten- 
tial. 

In Fig. @]we illustrate the time dependence of the edge 
flux for different barriers during a short time scale for 
0) = ifo(x), where the lines shift to the right by 
lowering the barrier. It can be seen from both panels 
of Fig. |4] that J(a\,t) and J(a,2,t), basically, increase 
during a short time interval and then reach at a charac- 
teristic time t p i a rough plateau, which shows a smoother 
behavior for J(a2, t), while J(ai, t) dives for a very short 
time to a negative value and oscillates more strongly. 
As a result, no feature can be clearly distinguished as 
a precise instant when tunneling begins. By using the 
continuity equations (|13[) . the decay time i; can be reli- 
ably estimated as 1/J p i(ci2), where J p i(a2) is the typical 
value of the flux at the plateau. The time scale when 
J(o2, t) develops a plateau also becomes larger, although 
it cannot be defined precisely. A crude estimate for the 
scale at which the plateau forms is t p \ ~ n/(2Eo), ap- 
proximately a factor of 2 less than the oscillation period 



period Tk determines the prefactor of the escape time as 
discussed below Eq. ([9]). Therefore, classically speaking, 
the decay by tunneling as a steady process begins when 
the electron hits the barrier. 

The inset in Fig. SJb) shows the behavior of the density 
at the edge. Contrary to J(ai, 0) = J(a,2, 0) = 0, p(a\, 0) 
and p(a2,0) are not zero, as the ground state for the 
potential U±(x) is not fully localized in the well, and the 
plateau in /9(<Z2,0) is clearly seen. 

Even though the numerically estimated times arc much 
shorter, a similar trend in variation with the height of 
the barrier is shown by the Biittikcr-Landauer time (BL 
time) , provided that nd remains unchanged. The traver- 
sal time of Biittiker and Landauer [lol ]. £bl = d/n, given 
by the barrier width d divided by the " semiclassical" ve- 
locity k = [2(Uo — Eo)] 1 / 2 , is an important time scale, 
especially in opaque conditions. With the conservation 
of nd, tsL is proportional to 1/k 2 . 

The dependence of the flux at the right edge on time 
with different width but keeping the height of the barrier 
is shown in Fig. 03a). It is interesting to find that the 
time that the flux begins to form a plateau is almost 
equal for the different widths, although the flux value 
at the plateau changes strongly, roughly as exp(— 2nd). 
In Fig. EJb), this characteristic time t p \ also remains 
almost unchanged for different Uq by keeping the same 
width ci2- We conclude that the observed scale of ir/ (2Eq) 
is universal and does not depend on the details of the 
potential. 

For comparison, we illustrate in Fig. 6 the flux at a 
short time scale for a transparent barrier nd = 0.25. As 
a matter of fact, the flux and the density have similar 
profiles. They grow from initial to the maximum value 
and decrease rapidly, contrary to the opaque behavior 
without an obvious peak. The peak is positioned at t sa 
0.4, similar to the time of plateau development in FigJS] 

For the initial state &(x, 0) = v / 2sin(7ra;) the flux 
is enhanced by an order of magnitude and oscillates 
stronger compared to the initial ^(x, 0) = ipo(x), be- 
cause sin(7rx) contains larger contributions from different 
cigenstatcs. Different from that with ^{x, 0) = <fo(x), 
the density is zero at t = outside the barrier. The 
time-dependence of the density at the edges with the 
initial state ^(x,0) = y/2sm(-Kx) is presented in Fig. 
[7] for a typical barrier. For various system parameters, 
the delay time between the maximum of the density at 
the right and the left edges is in a good agreement with 
the Biittiker-Landauer time £bl for these potentials. We 
conclude that £bl manifests itself as a delay time if over- 
the-barricr motion becomes essential due to the choice 
of the initial state. The momentum components that 
matter at first are the ones larger than y/2Uo, as the 
momentum distribution of sin(7rx) is broad. Moreover, 
its average local velocity at the right edge of the barrier 
UB(ci2,i) = J(a 2 ,t)/ p(a 2 ,t) decreases from a large value 
with some oscillations to a relatively stable one close to 
s/2Eq during a short interval, after which the real tun- 
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FIG. 6. (Color online) Time evolution of the flux at the right 
edge, where the solid line represents the flux, when Uo = 16 
and a,2 — 1.05. The dashed line corresponds to the flux for 
L^o = 10 and = 1.06. The two barriers are transparent as 
nd = 0.25, and the initial state is ipo{x). 



FIG. 7. Density evolution with the barriers Uo = 10, a? = 
1.63 (solid line at 02 and dashed line at ai). As the density 
at a,2 is relatively small, it is multiplied by the factor of ten. 
The barrier is opaque with nd = 2 and the initial state is 
*(z,0) = \/2sin(7rir). 



neling starts. This means that the forerunners just go 
above the barrier instead of tunneling. By contrast, the 
tunneling process for ^(x, 0) = <po(x) occurs indeed from 
the instant that the decay initiates, because VB(fl2,£) is 
always smaller than ^/2Uq. 

As discussed in [H, HI] with analytical models, tsL 
is a characteristic time describing different phenomena, 
among them over-the-barrier transients. This is rather 
paradoxical, considering its association with "tunnelling" 
in the defining formula, and has surely not been fully ap- 
preciated. A more intuitive understanding of this unex- 
pected role is still missing. 



C. Long-term dynamics 

We have now established that the short time dynam- 
ics of the probability flux is governed by robust time 
scales. We next investigate the long term dynamics with 
the goal to find out whether a suitable extrapolation of 
the long term scattering data can be used to gain infor- 
mation on the short time dynamics. We find that af- 
ter formation of the steady tunneling, the particle prob- 
ability density shows two distinct features. The first 
one is an almost uniform change in p(x, t) for x < a\ 
with p{x,t) ~ p{x, 0) exp(— t/ti). The second one can 
be viewed as a broad bump (wave packet) with relatively 
small density propagating with the velocity close to ^/2Eq 
and spreading in time. In this subsection, we consider the 
dynamics of the bump at long time scale and trace it to 
short times to obtain the operationally defined tunneling 
times. 

We assume that the flux J(X, t) and the density p(X 7 t) 
are measured by a detector located at X ^S> a 2 . It is 
shown in Fig. [5] for $>(x,0) — (fo(x), that the density at 
X 3> ai is nearly zero up to some time, then grows to 



a sharp maximum, and then decreases at timescales of 
ti with the sequential oscillations due to the diffraction 
in time phenomenon (2fj| . The profiles of the flux and 
density are very similar, with J(X,t) « y / 2E Q p(X, t). 

In the attosecond experiment analysis (3~H the tunnel- 
ing time was defined as the delay between the time when 
the barrier begins to be lowered and the time when elec- 
tron experiences the acceleration by external field, as can 
be extracted from the long-time behavior. Similarly, we 
can use the operational phenomcnological approach to 
define a tunneling time as 

*tun,i =tx~ ~ Q2 ) (15) 
v p 

where tx is the time when the detector measures the 

strongest first peak of the flux, v p = \J 2Eo is the velocity 
with which the particle moves out of the potential, and 
the energy 

/>oo 

Eo = E Q - / U (fi(x)dx (16) 

J a 2 

is slightly less than Eq , as the potential changes suddenly 
from Ui(x) to U2{x). As tx is a measurable quantity 
the tunneling time can be calculated by Eq. (fT5")) . For 
instance, for ^(x,t = 0) = (fo(x), Uq — 16, 0,2 = 1.4, 
when the detector is at X = 120, then tx = 50.3, and 
resulting ttun.i = 5.62. Another approach, similar to Eq. 
(fl~5)) . is to calculate the tunneling time as 



itun,2 =tx , (17) 

v x 

where vx is the velocity of the flux peak. For example, 
vx = 2.517 as defined by the motion from X = 120 to 
X = 122. According to Eq. (fTT)) . the corresponding 
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FIG. 8. Time evolution of density p(X, i) at given positions 
X = 60,90,120, provided that *O,0) = tp (x), U = 16, 
a-2 = 1.4. A strong initial peak is followed by oscillations due 
to the diffraction in time process. 



tunneling time is itun,2 = 3.18. As it can be assumed 
that the velocity of the electron is constant outside the 
barrier in a classical manner, we can conclude that the 
tunneling time extracted from the measurements by a 
remote detector is A-indcpcndcnt. Both Eq. (fT5|) and 
Eq. (fTT)) extrapolate electron motion from distance to 
the exit of the tunneling process. However, these two 
times are not equal either to the time of formation of the 
steady tunneling in Fig. |4]or to the decay time U. 

Fig. IH1 depicts how the flux evolves with time and dis- 
tance at two timescales for ^(x, 0) = (fo(x). The cur- 
vature of the maximum flux region in Fig. |9ja) shows 
that it takes some time for the flux to develop a constant 
speed in free space. In addition, Fig. Ela) demonstrates 
that a peak of the flux at a 2 appears at t « 1, smaller 
than ttun.2 m Eq. (fTT|) . 

It is expected in some models that the wave function 
at long times and distances can be obtained with the 
stationary phase approximation. Eq. ([2]) in this limit 
can be recast as 




(b) 45 50 55 

t 



FIG. 9. (Color online) The flux versus time and position 
after a short time near the barrier (a) and after a long time 
far away from the barrier (b). All the parameters are the same 
as in Fig. [8] The dashed straight lines serve as the guide for 
the eye only. The short-dash line in (a) demonstrates the 
nonlinearity of the peak position as a function of time. The 
long-dash line in (b) shows that this dependence is very close 
to linear at these times and distances. 



> a 2 ,t) 

1 



(18) 



2wi 



dk G(k) exp 



6(h) + kx 



fc 2 tl 
~2 



and the phase 9{k) + kx — k 2 t/2 can be expanded 
near the stationary point K satisfying the equation 
(d0(k) I ' dk)k=K+x— Kt = 0. However, G(k) in our calcu- 
lations is not a sufficiently smooth function near the K- 
points due to the resonances shown in Fig. [2j Therefore, 
the stationary phase approximation does not provide a 
satisfactory description of the peak propagation. 

Another factor that affects the tunneling time is the 
closing time ioj when the potential turns back from [^(x) 
to U\(x) and the tunneling is interrupted. The time evo- 
lution of the flux observed at the same remote position 
of detector with different closing times is demonstrated 



in Fig. 1 1 01 and ti is the exponential decay time intro- 
duced above. The peak of the flux with shorter closing 
time to appears earlier than that with long er to, as a 
result of the increased energy spreading j4Q| . After the 
closing time to, the probability to find the particle in the 
well remains almost unchanged. This shows that com- 
ponents with larger momenta tunnel through the barrier 
first, and therefore, closing the barrier can decrease the 
operationally defined times ttun,i and £tun,2- The mea- 
surement of the total number of particles escaped from 
the potential as the function of the closing time can help 
to find the timescale of the flux formation i p i: with the 
increase of to through this region the number of escaped 
particles which is quadratic in to for short times changes 
at to > t p i to a linear dependence. 
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t 



FIG. 10. (Color online) Flux as a function of time detected at 
X = 120 with different closing time to = oo (solid), to — ti/3 
(dashed), to — ti/6 (dot-dashed), to = ti/9 (dotted), ti is 
the life time for the particle in the case of to = oo. Other 
parameters are the same as those in Fig. [8] 



IV. CONCLUSIONS 

Motivated by recent experiments on the ultrafast atom 
ionization by optical fields, we have studied numerically 
exactly, for a rectangular barrier, the time-resolved tun- 
neling for short time scales and the further propagation 
at long times of an initially localized particle. The bar- 
rier we consider is opaque, however, not extremely so, 
to ensure a reasonable tunneling probability. The prob- 
ability density evolution on a short time scale much less 
than the decay time, depends strongly on the initial state. 
Dependent on how this state is prepared, this short-term 
motion can include both under-the-barrier tunneling and 
over-the-barrier propagation, as seen in the evolution of 
the density and flux at the barrier edges. The tunnel- 



ing starts instantly, however, some time of the order of 
tt/(2E ), where E is the ground state energy in the ini- 
tial potential, is required to develop the outgoing flux 
eventually proportional to the characteristic exponential 
decay rate 1/t;. As expected, there is a time delay be- 
tween the flux development at the left and the right edges 
of the barrier. If the initial state is the ground state of the 
potential at t < 0, the timescale of the flux development 
is much longer than the Buttiker-Landaucr traversal time 
expected for the given barrier parameters. However, if 
the initial state is more tightly localized, the Biittiker- 
Landauer time manifests itself as a time delay of the flux 
and density maxima between the left and right edges of 
the barrier. 

At long times we have considered the propagation of 
the escaped wave packet at distances much larger than 
the scale of the potential. From the operational definition 
of the escape time, related to the position of the maxi- 
mum of the wave packet density, we have estimated the 
time the particles escaped from the potential. This time 
is also much longer than the traversal time for the given 
barrier. To determine the effect of the time dependence 
of the barrier, we implemented escape time windows con- 
siderably shorter than the tunneling decay time t\. The 
increasing importance of faster components for shorter 
time windows leads to the extrapolated time estimated 
for the closing potentials smaller than that for the po- 
tentials permanently open. 
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